Combinatorial polynomials as moments, 

Hankel transforms and exponential Riordan 

arrays 



o 

(N 



Paul Barry 
School of Science 
Waterford Institute of Technology 
Ireland 
pbarryQwit . ie 



o 
u 



> 

o 

cn 

in 
o 



Abstract 

In the case of two combinatorial polynomials, we show that they can exhibited as 
moments of paramaterized families of orthogonal polynomials, and hence derive their 
Hankel transforms. Exponential Riordan arrays are the main vehicles used for this. 



1 Introduction 

Let [n] = 1,2, ... ,n, and let SP„ be the set of set-partitions of [n]. For a set-partition 
TT e SP„, let |7r| be the number of parts in vr. Then the n-th exponential polynomial, also 
known as the n-th Touchard polynomial (and sometimes called the n-th Bell polynomial 
[23]), is given by 



k=0 



where 



5(".*) = fE(-i)'-'(*u- 



j=0 

is the general element of the exponential Riordan array 

[l,e^'-l]. 

This is the matrix of Stirling numbers of the second kind A008277 , which begins 
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It is well known [11, 15, 18] that the Hankel transform of these polynomials is given by 

n 
k=l 

Now let 

k 

A{n,k) = Y,{-iy{k-j) 

be the general term of the triangle of Eulerian numbers. The matrix of these numbers 
A008292 begins 
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A[n, k) is the number of permutations in (3„ with k excedances. The Eulerian polynomials 
E\Jniz) are defined by 

n 

EUr,{z) = J2Mn,k)z\ 

k=0 

It is shown in [18] that the Hankel transform of these polynomials is given by 

n 
k=l 

These two results are consequences of the following two theorems. 

Theorem 1. The polynomials en{z) are moments of the family of orthogonal polynomials 
whose coefficient array is given by the inverse of the exponential Riordan array 

[e"(""-^\ - 1] . 

Theorem 2. The polynomials EU„(z) are moments of the family of orthogonal polynomials 
whose coefficient array is given by the inverse of the exponential Riordan array 

e^^'(l — z) — e^^ 

Note that in the case z = 1, the above matrix is taken to be \-r—^ t^I , whose inverse is 

' L l—x ' 1— X J ' 

the coefficient array of the Laguerre polynomials [2]. 

While partly expository in nature, this note assumes a certain familiarity with integer 
sequences, generating functions, orthogonal polynomials [5, 10, 22], Riordan arrays [17, 21], 
production matrices [9, 14], and the integer Hankel transform [4, 6, 13]. Many interesting 
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examples of sequences and Riordan arrays can be found in Neil Sloane's On-Line Encyclo- 
pedia of Integer Sequences (OEIS), [19, 20]. Sequences are frequently referred to by their 
OEIS number. For instance, the binomial matrix B ("Pascal's triangle") is A007318 . 
The plan of the paper is as follows: 

1. This Introduction 

2. Integer sequences, Hankel transforms, exponential Riordan arrays, orthogonal polyno- 
mials 

3. Proof of Theorem 1 

4. Proof of Theorem 2 



2 Integer sequences, Hankel transforms, exponential 
Riordan arrays, orthogonal polynomials 

In this section, we recall known results on integer sequences, Hankel transforms, exponential 
Riordan arrays and orthogonal polynomials that will be useful for the sequel. 

For an integer sequence a„, that is, an element of Z^, the power series fo{x) = J2T=o '^^^^ 
is called the ordinary generating function or g.f. of the sequence, while fe{x) = X^fc^o 1^^^ 
called the exponential generating function or e.g.f. of the sequence. a„ is thus the coefficient 
of in fo{x). We denote this by a„ = [x"']fo{x). Similarly, a„ = n![x"]/e(x). For instance, 
Fn = [x"] ^_^_^2 is the n-th Fibonacci number A000045 , while n\ = n\[x'^]j^, which says 
that is the e.g.f. of n\ A000142 . For a power series /(x) = Yl'^=o^"-^"^ with /(O) = 
and /'(O) 7^ we define the reversion or compositional inverse of / to be the power series 
f{x) = /t~^^(x) such that f{f{x)) = x. We sometimes write / = Rev/. 

The Hankel transform [13] of a given sequence A = {cq, ai, 02, •••} is the sequence of 
Hankel determinants {/iq, hi, /i2, • • • } where hn = |aj+j|"j^Q, i.e 



A — {an}n£No h — {hn}n 



ai a2 

O'n 0-n+l 



an 

0'2n 



The Hankel transform of a sequence a„ and its binomial transform are equal. 
In the case that a„ has g.f. g{x) expressible in the form 



ao 



1 — UqX 



1 — aix 



1 — a2X 



I32X^__ 

_Jhx^ 
1 — a^x - 
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(with (3i^0 for all i) then we have [11, 12, 24] 

n 

K = a^pr'pr' ■ ■ ■ Pi-iPn = n /^r'^'- (2) 

k=l 

Note that this is independent from an- In general and /3„ are not integers. Such a 
continued fraction is associated to a monic family of orthogonal polynomials which obey the 
three term recurrence 

Pn+l{x) = (X - an)Pn{x) - /3„p„_l(x), Pq{x) = 1, Pi{x) = X - Oq. 

The terms appearing in the first column of the inverse of the coefficient array of these 
polynomials are the moments of family. 

The exponential Riordan group [1, 7, 9], is a set of infinite lower-triangular integer matri- 
ces, where each matrix is defined by a pair of generating functions g{x) = go + gix + g2x'^ + . . . 
and f{x) = fix + /2X^ + . . . where /i ^ 0. The associated matrix is the matrix whose i-th 
column has exponential generating function g{x)f{xy/i\ (the first column being indexed by 
0). The matrix corresponding to the pair /, g is denoted by [g, /]. It is monic if g^ = 1. The 
group law is given by 

[gJ]*[h,l] = [g{hof),lof]. 

The identity for this law is / = [1, x] and the inverse of [g, f] is [g, f]~^ = [1/(5' o /), /] where 
/ is the compositional inverse of /. We use the notation eTZ to denote this group. If M is the 
matrix [g, /], and u = {un)n>o is an integer sequence with exponential generating function U 
(x), then the sequence Mu has exponential generating function g{x)U{f{x)). Thus the row 
sums of the array [g, f] are given by g{x)e^^^'> since the sequence 1,1,1,... has exponential 
generating function e^. 

Example 3. The binomial matrix is the matrix with general term (^). It is realized by 
Pascal's triangle. As an exponential Riordan array, it is given by [e^,x]. We further have 

([e^x]^ = [e'"^x]. 

Example 4. We have 

[e-(e^-i),e"-l] = [e'^'^~^\x] ■[l,e"-l]. 
A more interesting factorization is given by 
Proposition 5. The general term of the matrix L = \e^'^'^'^ — l] is given by 

i^n,. = X^5(r^,j)(^)^'-'- 

Proof. A straight-forward calculation shows that 

[e-(<=^-i),e"-l] = [l,e"-l] ■ [e^",x]. 

The assertion now follows since the general term of [1, — 1] is S'(n, k) and that of [e^^, x\ 
is 0^""^ □ 
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As an example of the calculation of an inverse, we have the following proposition. 
Proposition 6. 



[e-(-^-i),e"-l] ' = [e-"",ln(l + 
Proof. This follows since with 



X 



we have 



Proposition 7. 



fix) = e'-l 
f{x) = ln(l + x) 



□ 



e^^fl — z) — e^^ 



1 -1 



e^^ — ze^ e^^ — ze 
Note that in the case z = 1, we have 

1 



1 + zx, 



z-1 



In 



1 + zx 
1 + x 



X 



1 — X 1 — X 



X 



1+x 1+x 



Proof. This follows since with 



we have 



e^^d - z) 



z-l 



In 



1 + zx 
1 + x 



□ 



An important concept for the sequel is that of production matrix. The concept of a 
production matrix [8, 9] is a general one, but for this note we find it convenient to review it 
in the context of Riordan arrays. Thus let P be an infinite matrix (most often it will have 
integer entries). Letting ro be the row vector 

ro = (1,0,0,0,...), 

we define r^ = rj_iP, i > 1. Stacking these rows leads to another infinite matrix which we 
denote by Ap. Then P is said to be the production matrix for Ap. If we let 



M^ = (1,0,0,0,...,0,...) 



then we have 



At 



u^P 
u^P^ 

\ J 



and 

DAp = ApP 

where D = (^ij+i)ij>o (where S is the usual Kronecker symbol). In [14] P is called the 
Stieltjes matrix associated to Ap. In [9], we find the following result concerning matrices 
that are production matrices for exponential Riordan arrays. 

Proposition 8. Let A = (an,fc)„ fc>o = [fi'(^))/(^)] exponential Riordan array and let 

c{y) = Co + ciy + C2?/^ + . . . , r{y) = ro + riy + r2y^ + . . . (3) 
he two formal power series that that 

r{f{x)) = fix) (4) 

cifix)) = (5) 
9{x) 



Then 



^^'Cjfln.i (6) 



(ii) an+i,k = roan,k-i + j;^^iKci-k + kri_k+i)an,i (7) 



i>k 



or, defining c_i = 0, 



0'n+l,k 



i>k-l 



Conversely, starting from the sequences defined by (3), the infinite array {,0'n,k)nk>o defined 
by (8) is an exponential Riordan array. 

A consequence of this proposition is that P = {pi,j)ij>o where 

Pi J = -^{Ci-j + JTr-j+l) (C_l = 0). 

Furthermore, the bivariate exponential generating function 

(ppit.z) = ^Pn,k't^ 
n,k 

of the matrix P is given by 

0p(t,z) = e''{c{z) +tr{z)). 

Note in particular that we have 

r{x) = f{f{x)) 

and 

, _ g'ifjx)) 



n! 



6 



Example 9. We consider the exponential Riordan array A094587 . This array [2] 

has elements 



1 
1 
2 
6 

24 
120 




1 
2 
6 

24 
120 





1 
3 

12 
60 






1 
4 
20 











1 
5 1 



V 



and general term [k < n]y with inverse 



/ 1 

-1 





V '■■ 





1 

-2 









1 

-3 









1 









1 
-5 1 



which is the array [1 — x,x]. In particular, we note that the row sums of the inverse, 
which begin 1, 0, —1, —2, —3, . . . (that is, 1 — n), have e.g.f. (1 — x) exp(x). This sequence 
is thus the binomial transform of the sequence with e.g.f. (1 — x) (which is the sequence 
starting 1, —1, 0, 0, 0, . . .). In order to calculate the production matrix P of [j^, x] we note 
that /(x) = X, and hence we have f'{x) = 1 so f'{f{x)) = 1. Also g{x) = leads to 
g'{^) = (i-x)i ) so, since f{x) = x, we get 

Thus the generating function for P is 



1 



X 



1-z 



Thus P is the matrix [yz;;,^^] with its first row removed. 

Example 10. We consider the exponential Riordan array [1, ^3^]. The general term of this 
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matrix [2] may be calculated as follows: 



T 

J- 7 



n.k 



kV 



n 



X] 



K ! 



\x 



n—k' 



j=0 

oo 



k+j 



x-" 



n 
k\ 

3=0 

n\fk + n — k — 1 
k\\ n — k 
n\ fn — 1 
k\ \n — k^ 

Thus its row sums, which have e.g.f. exp (yz^); have general term X^fc=o f[ (n-D ■ This is 
A000262 , the 'number of "sets of lists" : the number of partitions of {1, .., n} into any number 
of lists, where a list means an ordered subset'. Its general term is equal to (ra — l)!L„_i(l, — 1). 
The inverse of [l, ^3;] is the exponential Riordan array [l, ? A111596 . The row sums 
of this sequence have e.g.f. exp (3^), and start 1, 1, —1, 1, 1, —19, 151, . . .. This is A111884 . 
To calculate the production matrix of [l, we note that g'{x) = 0, while f(x) = with 
fix) = Thus 



x) 



f'ifi^)) = (1 , , 
and so the generating function of the production matrix is given by 



The production matrix of the inverse begins 



/ 1 



2 1 

2 4 

6 









1 
6 

12 









8 1 
20 10 






1 
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Example 11. The exponential Riordan array A 



( 1 
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1 
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1 











6 


18 


9 


1 








24 


96 


72 


16 


1 





120 


600 


600 


200 


25 


1 



\ 



has general term 



T 



n\ In 



Its inverse is j+^j with general term (-1)" ''flQ)- This is A021009 , the triangle of 
coefficients of the Laguerre polynomials Ln{x). The production matrix [j^zjyz"] is given 
by 



/I 1 

1 3 

4 

9 











1 

7 

16 








1 
9 

25 








1 

11 



y : : : : : : •• / 
Example 12. The exponential Riordan array [e^,ln (j^)], or 
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.. \ 
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1 
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1 








. 
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1 





. 
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24 


29 


10 


1 


. 




1 


89 


145 


75 


15 


1 . 



















is the coefficient array for the polynomials 



-n, X] —1] 



which are an unsigned version of the Charlier polynomials (of order 0) [10, 16, 22]. This is 
A094816 . It is equal to 

[e^,x] l,ln 



X 
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or the product of the binomial array B and the array of (unsigned) Stirling numbers of the 
first kind. The production matrix of the inverse of this matrix is given by 



/-II 
1-21000 

2-3100 

3 -4 1 

4 -5 1 

5 -6 

V ; ; ; ; ; ; 



which indicates the orthogonal nature of these polynomials. We can prove this as follows. 
We have 

1 



e^,ln 



1 — X 



Hence g{x) = e ^ and f{x) = 1 — e ^ . We are thus led to the equations 

r(l - e"^) = e~^, 
c(l-e-") = -e"", 

with solutions r(x) = 1 — x, c{x) = x — 1. Thus the bivariate generating function for the 
production matrix of the inverse array is 



which is what is required. 



3 Proof of Theorem 1 

Proof. We show first that with L = \e^'^'^''~^\ — l] , the matrix L^^ which is given by 

L^i = [e"(^'-i), - 1] = [e-"^ ln(l + x)] , 

is the coefficient array of a family of orthogonal polynomials. To this end, we calculate the 
production array of [e^^^''~^\ — l] . We have f{x) = e^ — 1, f'{x) = and f{x) = In(H-x). 
Thus 

c{x) = f'{f{x)) = l + x. 
Similarly, for g(x) = e^^^''~^\ we have g'{x) = ze^^^'°~^'^~^^ and so 

g'ifjx)) ze^^{l + x) 
r(x) = = = z(l + x). 

Thus the production matrix sought has generating function 

e'^^iciw) + tr{w)) = e*"'(l + w + t{z{l + w))). 
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Thus the production array Pl is tri-diagonal, beginning 



/ 


z 


1 


















z 


^ + 1 


1 


















22 


2 + 2 


1 


















32 


2 + 3 


1 


















^z 


2; + 4 


1 


















52 


2 + 5 . 




V 

















Now it is well known that 



E 

/c=0 



fc 2(e^-l) 



and hence the polynomials e„(2;) are the moments the family of orthogonal polynomials 
whose coefficient array is L~^. □ 

Corollary 13. The Hankel transform of en{z) is z^ ^ ^ 11^=1 

Proof. From the above, we have that the generating function of en{z) is given by the con- 
tinued fraction 

1 



1 — 2;x — 



zx 



1 - (2 + l)a; 



22X^ 



1- (2 + 2)x 



320;^ 
1 - ■• 



In other words, = nz. Thus the Hankel transform of en{z) is given by 



n n 

jn—k+1 

k=l k=l 



n-k+l _ J"+^) 



k=l 



□ 



We note that the Hankel transform of the row sums of L = ^e^^'^'' ^\e^ — l] is equal to 



it=i 



Note also that if we take z = e*, we obtain a solution to the restricted Toda chain [3]. 



4 Proof of Theorem 2 



Proof. We show first that with L 
L 1 = 



e^'=(l-z) 



the matrix L ^ which is given by 



X ' ^zx 



ze^ e 



ze^ 



^ , 1 / 1 + 2X 

1 + ZX, m 



2 — 1 V 1 + X 
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is the coefficient array of a family of orthogonal polynomials. To this end, we calculate the 



production array of L 
and 



e^^(l-z) 



fix) 



. We have /(x) 

1 - z)2gx(l+.) 



L,/(x) = ^ln(M|l) 



x\2 



Thus 
Also g{x) 



C(X) = /'(/») = (1+X)(1 + 2X). 

which implies that g\x) = xf'{x) and so 



r{x) 



z(l +x) 



Thus the generating function of Pl is given by 

e^'^iciw) + tr{w)) = e*'"((l + w){l + zw) + t{z{l + w))). 
Thus the production array Pl is tri-diagonal, beginning 





z 


1 


















z 


2z + l 


1 


















Az 


3^ + 2 


1 


















9z 


4z + 3 


1 


















IQz 


5z + 4 


1 


















2bz 


6^ + 5 . 




v 

















Now it is known that 



k=0 



k\ 



X 



e — ze 



and hence the polynomials en{z) are the moments the family of orthogonal polynomials 
whose coefficient array is L~^. □ 

Corollary 14. The Hankel transform of EU„(2;) is zi '^^ YYk=i 

Proof. From the above, we have that the generating function of EU„(2;) is given by the 
continued fraction 



zx- 



zx 



{2z + l)x 



Azx 



(32 + 2)x 



9zx 



In other words, /3„ = n'^z. Thus the Hankel transform of EU„(z) is given by 



2\n-k+l 



fc=i 



k=l 



k=l 



□ 
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